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A new combined cooling, heating and power (CCHP) system is proposed. This system is driven by solar 
energy, which is different from the current CCHP systems with gas turbine or engine as prime movers. 
This system combines a Rankine cycle and an ejector refrigeration cycle, which could produce cooling 
output, heating output and power output simultaneously. The effects of hour angle and the slope angle of 
the aperture plane for the solar collectors on the system performance are examined. Parametric opti¬ 
mization is conducted by means of genetic algorithm (GA) to find the maximum exergy efficiency. It is 
shown that the optimal slope angle of the aperture plane for the solar collectors is 60° at W a.m. on June 
12, and the CCHP system can reach its optimal performance with the slope angle of 45° for the aperture 
plane at midday. It is also shown that the system can reach the maximum exergy efficiency of 60.33% 
under the conditions of the optimal slope angle and hour angle. 

© 2009 Elsevier Ltd. All rights reserved. 


1. Introduction 

The trigeneration concept refers to the possibility of efficiently 
and profitably combining cooling, heating and power production 
(CCHP) to satisfy the consumption needs during the energy system 
operation. In a competitive energy market framework, the adoption 
of CCHP systems may become profitable with respect to traditional 
systems, where electricity, heat, and cooling are produced or 
purchased separately. 

There are several current CCHP technologies to perform the 
trigeneration according to the prime movers, such as steam turbine, 
reciprocating internal combustion engine, gas turbine, micro gas 
turbine, Stirling engine, fuel cell. Wang et al. [1] had made an 
extensive and an intensive review of CCHP in the literature. Wang 
et al. [2,3] also investigated the CCHP system driven by Stirling 
engine and gas turbine. Li et al. [4] conducted the energy utilization 
evaluation for the CCHP systems mainly with gas engine and gas 
turbine as prime movers. 

Although steam turbine, reciprocating internal combustion 
engine and gas turbine that can be considered as the conventional 
prime movers still make up most of the gross capacity being installed, 
and micro gas turbine, Stirling engine and fuel cell present a prom¬ 
ising future for prime movers in CCHP system, the current CCHP 
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technologies must inevitably consume the fossil fuels, which would 
disappear with energy consumption and economic development in 
the future. So, it is necessary to utilize renewable energy. 

The main objective of the present study is to propose a new 
combined cooling, heating and power (CCHP) system driven by 
solar energy. This system combines Rankine cycle and ejector 
refrigeration cycle, and produces cooling output, heating output 
and power output simultaneously. The ejector refrigeration cycle, 
which many studies [5-11] have been devoted to, is different from 
the absorption refrigeration cycle in most of CCHP systems. It has 
some advantages such as fewer movable parts and low operating, 
installation and maintenance cost. In addition, the ejector refrig¬ 
eration cycle has the possibility of using a wide range of refrigerants 
with the system. 

In the present study, the simulation of the CCHP system driven 
by solar energy is achieved using a simulation program. The effects 
of hour angle and the slope angle of the aperture plane for the solar 
collectors on the system performance are examined. In addition, 
the parameter optimization is achieved with exergy efficiency as 
the objective function by means of genetic algorithm under the 
given condition. 

2. System description 

The proposed CCHP system is driven by solar energy, and 
combines the Rankine cycle and the ejector refrigeration cycle, 
which could produce cooling, heating and power simultaneously. 
Fig. 1 illustrates the CCHP system. The overall system is divided into 
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Nomenclature 

Subscripts 



a 

acceptance half angle 

A 

surface area of solar collectors, m 2 

b 

beam 

b 

width of absorber, m 

br 

branching 

c 

heat capacity, kj kg -1 

B 

boiler 

C 

concentration ratio 

C 

condenser 

D 

diameter, m 

d 

diffuse radiation; diffuser section of ejector 

E 

exergy, kW 

e 

effective 

h 

enthalpy, kj kg' 1 

exg 

exergy 

H 

height, m 

extr 

extraction 

I 

hourly radiation, W m 2 

E 

evaporator 

k 

heat transferring coefficient, W m“ 2 K _1 

fi 

inlet of solar collectors 

L 

length of tube, m 

fo 

outlet of solar collectors 

m 

mass flow rate, kg s _1 

g 

global 

ma 

average number of reflections 

he 

average temperature of heater 

n 

the day of the year 

H 

heater 

N 

number of tubes 

i 

Inlet 

P 

pressure, MPa 

inst 

Instantaneous 

Q 

heat load, kW 

1 

uniform 

R 

tilt factor; ratio 

lo 

Loss 

s 

entropy, kj kg -1 K _1 

m 

mixing section of ejector 

S 

total effective flux, W m~ 2 

mf 

mixed fluid 

t 

temperature, °C; time 

n 

nozzle 

T 

temperature, K 

nl 

inlet of nozzle 

u 

velocity, m s _1 

n2 

outlet of nozzle 

U 

overall loss coefficient, W m -2 1< -1 

NET 

Net 

W 

width, m 

P 

pressure 



pf 

primary flow 

Greek letters 

PI 

pump 1 

a 

absorptivity of the absorber surface 

P2 

pump 2 

P 

slope angle of the aperture plane for solar collectors, ° 

s 

isentropic process 

5 

declination 

sc 

average temperature of solar collectors 

V 

efficiency 

sf 

secondary flow 

6 

angle of incidence 

thm 

thermal 

fi 

entrainment ratio 

T 

turbine 

P 

reflectivity of the concentrator surface; density, kg nrT 3 

u 

useful 

T 

transmissivity of the cover 

z 

zenith angle 

(0 

hour angle, ° 

0 

environment 

(P 

latitude, ° 




two subsystems: the solar collector subsystem and the CCHP 
subsystem. 


The hourly global radiation J g reaching a horizontal surface on 
the earth is given by: 

/ g = / b +/d (i) 


2.1. The solar collector subsystem 


where 


This subsystem consists of solar collectors, a thermal storage 
tank and an auxiliary heater. The solar collectors are used as 
a main energy source to supply heat to the system. The thermal 
storage tank is used as the thermal source when solar radiation is 
not sufficient. The auxiliary heater is installed as the backup 
energy source to boosts the temperature of thermal storage tank 
to the allowable reference temperature when the temperature of 
thermal storage tank drops below the allowable reference 
temperature. 

A compound parabolic collector (CPC) is used to collect the solar 
radiation because it achieves higher concentration for large 
acceptance angle and requires only intermittent sun-tracking. In 
addition, the CPC can achieve the higher temperature than flat plate 
collector for power generation. In the present study, the symmetric 
CPC with flat absorber is placed on a horizontal east-west axis and 
oriented with its aperture plane sloping at an angle of 45°. 


4 = 4nCOS 0 Z (2) 

where Jb n is beam radiation in the direction of the rays, and 0 Z is 
angle of incidence on a horizontal surface, i.e. the zenith angle. 

It is postulated that for a clear cloudless day: 


7 bn = A exp 


B 

cos 0 Z 


and 



4 = C/ bn (4) 

where A, B and C are constants [12]. 

Since CPC for absorbing radiation is tilted at an angle to the 
horizontal, it is necessary to calculate the tilt factor for radiation. 
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The ratio of the beam radiation flux falling on a tilted surface to 
that falling on a horizontal surface is called the tilt factor R b for 
beam radiation, which can be given by: 


_ cos 0 _ sin dsin (0 - /?) + cos <Scos wcos(0 - 0) 
b — cos 0 Z ~ sin <5sin 4> + cos <Scos wcos 4> 

where the declination 5 can be given by Cooper [13]. 



5 


23.45sin 


'360 

365 


(284 + n) 



The tilt factor for diffuse radiation is the ratio of the diffuse 
radiation flux falling on the tilted surface to that falling on a hori¬ 
zontal surface. The value of this tilt factor depends upon the sky 
dome seen by the tilted surface. Assuming that the sky is an 
isotropic source of diffuse radiation, we have 


ma 


1 

2sin 0 a 



(1 - sin 0 a )(l + 2sin 0 a ) 
2sin 2 0 a 


( 11 ) 


Thus, the effective reflectivity of the concentrator surface is 
given by: 



( 12 ) 


where p e is the effective reflectivity of the concentrator surface for 
all radiation, and p is the reflectivity value for a single reflection. 

Because of its large acceptance angle, a CPC accepts both beam 
and diffuse radiation. The beam radiation flux falling on the aper¬ 
ture plane is /t>Rt> while the diffuse radiation flux within the 
acceptance angle is given by (/d/C). Thus, the total effective flux 
absorbed at the absorber surface is: 


S 


4> R b + / 


Tp e a 


(13) 


fid - (7) 

since (1 +cos/?)/2 is the radiation shape factor for a tilted surface 
with respect to the sky. 

The geometry of a two-dimensional CPC is shown in Fig. 2. The 
concentration ratio is given by: 


W _ 1 

b ~ sin 0 a 



where 0 a is the acceptance half angle. 

The height to aperture ratio of the concentrator is given by: 


// = /fl + ■ C )cos (9) 

W 2 \ sin 6 a J 

Rabl [14] has shown that the ratio of the surface area of the 
concentrator to the area of the aperture is given by the expression: 


A 

WL 


sin0 a (l +sin0 a ) 


cos0 a 

sin 2 0 a 


-bin 


f (l+Sin(? a )(l+COS(?a) 

V2cos0 a 

lsin0 a 

cos0 a + (2 + 2sin0 a )°' 5 


(l+sintfa) 1 ' 5 . 


( 10 ) 


Rabl has also shown that the average number of reflections ma 
undergone by all radiation falling within the acceptance angle 
before reaching the absorber surface is given by the expression 


Transparent cover 
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The useful heat gain rate can be given by: 


Qu = FrWL 


U] 


lo 


c 


(Tb - To) 


where 



jnCpT 

bU l0 L\ 


and 


exp 


F'bU\ 0 L ] \ 

mc p J j 


1 

F 



1 jb 


The instantaneous collection efficiency is given by: 


(14) 


(15) 


(16) 


4inst 


Qu 

(/ 5 J ? 5 +WWt 


(17) 


The insulated liquid storage tank with sensible heat storage is 
shown in Fig. 3, receiving energy from an array of collectors and dis¬ 
charging energy to a load for use in the CCHP subsystem. We assume 
that the liquid in the tank is always well-mixed and consequently is at 
T\. T\ is defined as the water temperature in the tank after two streams from 
solar collector and boiler are well-mixed. An energy balance on the 
thermal storage tank yields the following equation [12]: 

[(pVc p \+(pVc p ) t }^ = Qu - Qi oad - (LM) t (T, - T 0 ) (18) 

where (pVc p )j represents the heat capacity of the liquid in the tank, 
(pVc p ) t , the heat capacity of the tank material, Qu, the rate of useful 
heat gain received from the collectors, Qj oad , the rate at which 
energy is being discharged to the load, (UA\, the product of the 
overall heat-transfer coefficient and surface area of the tank. The 
heat capacity term (pVc p ) t is likely to be of importance for small¬ 
sized tanks only. For large tanks, its value may be negligible in 
comparison to (pVc p ) { . Denoting the sum of the two heat capacities 
by the symbol (pVc p ) e and integrating differential Eq. (18) under 
the assumption that Qu,Qioad an d To are constants, we get, subject to 
the initial condition t = 0, 7j = Ti [12]: 


Qu - Qioad - mm - To) r m\t - 

Qu - Qioad - mm - T 0 ) (pvc p ) e . 


(19) 


From collectors t 0 i oad 



The values of Qi oa d may be calculated from the inlet and outlet 
temperature and flow rate as: 

Qioad = ^loadCp(Tl — ^i) (20) 


2.2. The CCHP subsystem 

The CCHP subsystem consists of a few major components: boiler, 
turbine, evaporator, heater, regenerator, condenser, ejector, pumps 
and throttle valve. The boiler is a device in which high pressure and 
temperature vapor is generated by absorbing heat from thermal 
storage system. The high pressure and temperature vapor are 
expanded through the turbine to generate power. The extracted 
vapor from the turbine is divided into two streams. One stream 
enters the heater to supply heat load to heat user. The other stream as 
the primary vapor enters the supersonic nozzle of ejector. The very 
high velocity vapor at the exit of the nozzle produces a high vacuum 
at the inlet of the mixing chamber and entrains secondary vapor into 
the chamber from the evaporator. The two streams are mixed in the 
mixing chamber. Then the mixed stream becomes a transient 
supersonic stream. On entering the constant cross-section zone, 
a normal shock wave occurs, accompanied by a significant pressure 
rise. After the shock, the velocity of the mixed stream becomes 
subsonic and decelerates in the diffuser. The stream from ejector and 
the turbine exhaust are mixed together, and they enter the 
condenser where they condense from vapor to a liquid by rejecting 
heat to the surroundings. One part of the working fluid leaving the 
condenser enters the evaporator after passing through the throttle 
valve, and the other part flows to the regenerator after passing 
through one pump. The stream with relatively high temperature 
leaving the heater heats the stream from a pump in the regenerator 
where the two streams are mixed together. The pressure of the liquid 
leaving the regenerator is increased to the boiler pressure by the 
other pump, and it is pumped to the boiler to be vaporized again. The 
working fluid entering the evaporator is at low pressure and 
temperature. In the evaporator, it is vaporized by absorbing heat 
from the cooled media. Thus, a cooling effect is produced. 

The ejector is the key component in this CCHP system. In the 
present study, the ejector performance simulation is carried out 
based on the one-dimensional constant pressure flow model used 
by most researchers in studying ejector refrigeration systems. The 
basic principle of the model was introduced by Keenan et al. [15] 
based on gas dynamics, and developed by Munday et al. [16], Huang 
et al. [17] and Ouzzane et al. [18]. The following assumptions are 
made for the analysis. 

(1) The flow inside the ejector is in steady state and one¬ 
dimensional. 

(2) The primary and secondary flows at the inlet of the ejector and 
the mixed flow at the outlet of the ejector are at stagnation 
conditions. 

(3) For simplicity, the effects of frictional and mixing losses in the 
nozzle, diffuser and mixing section are taken into account by 
using the nozzle efficiency, the mixing efficiency and the 
diffuser efficiency. 

(4) Whenever the mixed flow is supersonic, a normal shock wave 
is assumed to occur upstream of diffuser inlet. 

(5) Mixing process in the mixing section of ejector occurs at 
constant pressure and complies with the conservation of 
energy and momentum. 

(6) The ejector does not exchange heat with the surroundings. 

Based on the above assumptions, the mass, momentum and energy 
equations are applied to sections for the ejector and then the 
following equations are obtained. 
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In the nozzle section, the energy conservation equation for the 
adiabatic and steady primary flow is given as: 


frpf,n2 + U pf,n2/2 — ^pf,nl + U pf,nl 

Compared with the outlet velocity of primary flow u P f >n 2 , the 
inlet velocity of primary flow u P f >n i could be negligible. According to 
Eq. (21), the outlet velocity of primary flow u p f in 2 is expressed as: 


t*pf,n2 — \j^ (frpf,nl W) 


The nozzle efficiency is given as: 


( 22 ) 


Vn 


fopf,nl ~ ftpf,n2 
^pf,nl — ^pf,n2,s 


(23) 


hence 


t*pf,n2 — (^pf,nl ^pf,n2,s) 


(24) 




The entrainment ratio of ejector is given as: 


™sf 


™ P f 


(25) 


In the mixing section, the momentum conservation equation is: 


tftpft*pf ? n2 "T tn sf U s f n2 — (m p f + 
hence 


(26) 


where 


^pf,n2 + M*^sf,n2 

1 + /X 


(27) 


^sf,n2 — (/l s f ^ s f n2 ) 


(28) 


Compared with the primary flow velocity zi p f >n 2 , the velocity of 
secondary flow zj s f >n 2 is negligible. Then Eq. (27) becomes: 


_ u pf,n2 
mf.m.s - -[ + ^ 

The mixing efficiency is given as: 


(29) 


u 


Vm 


mf.m 


U 


(30) 


mf.m.s 


The average velocity of mixed flow is written according to Eqs. 
(29) and (30) as: 


t*mf,m — t/pf n2 V^?m/(1 T- /x) 


(31) 


The energy conservation equation for the mixing section is: 


m pf V,n2 + 


u 


pf,n2 


+ m s f I h s f n2 + 


u 


sf,n2 


(rripf + m S f) f h m f,m T- 


u 


mf.m 


(32) 


The enthalpy of mixed flow is written according to Eqs. (22), 
(28) and (32) as: 


(21) ^mf,m — (^pf,nl + M^sf)/(1 + M) ^ m f m /2 


(33) 


In diffuser section, the mixed fluid converts the kinetic energy into 
pressure energy. Assuming the compression process is isentropic, 
the energy equation for the diffuser section is: 


2 (^mf,m ^mf,d,s) ^mf,d,s ^mf, 


m 


(34) 


The actual energy equation for the diffuser section is: 


~(u 2 

2 y u mf,m 


U 


mf,d 


— ^mf,d ^mf.m 


(35) 


The diffuser efficiency is given as: 


Vd 


^mf,d,s ^mf.m 
^mf.d — ^mf.m 


(36) 


Compared with the velocity u m f >m , the outlet velocity of the mixed 
fluid u m f t d can be neglected, the actual outlet enthalpy of the mixed 
fluid /i m fd is written according to energy Eq. (35) as: 


^mf,d — ^mf,m ^mf.m/^ 


(37) 


In addition, according to Eq. (36), the actual outlet enthalpy of 
the mixed fluid is also expressed as: 


^mf,d — ^mf.m + (^mf,d,s ^mf.m) / Vd 


(38) 


It is known that the performance of ejector is evaluated by its 
entrainment ratio / x . Based on above Eqs. (24), (31) and (38), it can 
be derived as: 

M — \JVnVmVd (^pf,nl — ^pf,n2,s) j (^mf,d,s — ^mf,m) — 1 (39) 

When the inlet state parameters of primary flow, secondary flow 
and back pressure of the ejector are given, the value of entrainment 
ratio could be found using iterative calculation. 

Extraction mass flow ratio is defined as the extraction mass flow 
divided by the turbine inlet mass flow, given by: 


^extr — 


Hlextr 

rrij 


(40) 


Mass flow branching ratio is defined as the primary mass flow of 
ejector divided by the extraction mass flow, given by: 


R 


m 


br 


pf 


t^extr 


(41) 


The basic equations obtained from the conservation law for 
energy in the components are written as follows: 


For evaporator : Qe = m E (h 9 - h 8 ) 
For boiler : Qb = m B (h 2 - h\) 

For condenser : Qc = m c (h 14 - h 5 ) 


(42) 

(43) 

(44) 


Forturbine: Wj = mj(h 2 -h 3 ) + (mj-m extr )(h 3 -h-i 3 ) (45) 
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For pump 1 : W P1 = m P i(hn - h 6 ) (46) 

For pump 2 : W P2 = m P2 (/ii - h 12 ) (47) 

For heater : Qh = m H (/i 3 - h 10 ) (48) 

For regenerator : m P2 h 12 = m H /i 10 + m P1 hn (49) 


2.3. System performance 


The overall system performance can be evaluated by the 
thermal efficiency and exergy efficiency. Thermal efficiency of the 
overall system is defined as the useful energy output divided by 
the total solar energy input, given by: 


^thm 


Wnet + Qe + Qh 

Qu 


(50) 


where W NET is the power output from the turbine, reduced by the 
power input to the pump, Qe is the cooling output, Qh is the heat 
output and Qu is the total heat added to the system from the solar 
collectors. 

From the viewpoint of the first law of thermodynamics and 
energy conservation used to determine the overall thermal effi¬ 
ciency, work and heat are equivalent. On the other hand, exergy, 
based on the second law of thermodynamics, quantifies the 
difference among work, heat and cooling in terms of irreversibility, 
or change in energy quality. Therefore, the exergy efficiency is 
chosen to be the criterion for the overall system performance 
evaluation. 

Exergy efficiency is defined as the exergy output divided by the 
exergy input to the overall system, which can be given by: 


Vexg 


W N et + Ee + £h 



(51) 


E u is the exergy input to the overall system from the solar 
collector, which is given as: 


Eu = Qu (l -(52) 

E e is the exergy associated with the refrigeration output, which is 
calculated as the working fluid exergy difference across the 
evaporator. 


Ee = m E' [(^E,in ~ ^E,out) ~ ^o( 5 E,in “ s E,out)_ 

Eh is the heat exergy output, which is given by: 

£ " -«"('-£) 


(53) 


(54) 


3. Methodology 


is 0.89, reflectivity of concentrator is 0.87, and absorptivity of 
absorber surface is 0.94. The number of CPC arrays is assumed to be 
1000, that can be designed according to the requirement load for the 
different applications. 

Water is used as storage liquid under the pressure of 0.5 MPa in 
the sensible heat storage system. The cylindrical hot water storage 
tank, 0.8 m in diameter and 1.0 m high, is made from a steel plate 
(p = 7800 kg m -3 , c p = 0.46 kj kg -1 K _1 )0.006 m thick. And it is 
insulated all round with glass wool insulation 0.02 m thick. 

R123 is selected as the working fluid because it is recognizable 
as a low pressure refrigerant that is nontoxic, nonflammable and 
non-corrosive. 

For the CCHP system simulation, the following assumptions are 
also made: 

(1) The system reaches a steady state, and pressure drop in pipes 
and heat losses to the environment in the boiler, turbine, 
heater, ejector, condenser, regenerator and evaporator are 
neglected. 

(2) The flow across the throttle valve is isenthalpic. 

(3) The condenser outlet state is saturated liquid, and its temper¬ 
ature is assumed to be approximately 5 °C higher than the 
environment temperature. 

(4) The working fluid at the evaporator outlet is saturated vapor. 

The main assumptions for the CCHP system simulation are 
summarized in Table 1. 

The efficiencies of the ejector are assumed to be p n = 0.90, 
rj m = 0.85, pa = 0.85 by making reference to literatures [8,17,19]. 

In order to improve the performance of the CCHP system, it is 
necessary to conduct a parametric optimization to reach the best 
performance. In the present study, the exergy efficiency, which can 
evaluate the performance of the CCHP system, is selected as 
objective function for parameter optimization. 

Parameter optimization is achieved by means of genetic algo¬ 
rithm to reach the maximum exergy efficiency. The GA, which is 
presented firstly by professor Hooland [20], is a stochastic global 
search method that simulates natural biological evolution. Based on 
the Darwinian survival-of-fittest principle, the genetic algorithm 
operates on a population of potential solutions to produce better 
and better approximations to the optimal solution. The GA differs 
from more traditional optimization techniques because it involves 
a search from a population of solutions and not from a single point. 

The GA encodes a potential solution to a specific domain problem 
on a simple chromosome-like data structure (which constitutes an 
individual), where genes are parameters of the problem to be 
solved. In the present study, the float-point coding is used in 
parameter optimization for the CCHP system. Each chromosome 
vector is coded as a vector of floating point numbers of the same 
length as the dimension of the search space. Chromosome is defined 
as a real number vector, X = (xi, x 2 ,..., x n ), x^R, i = 1,2,..., n, where i 
is i-th parameter for the CCHP system, n is the number of optimizing 
parameters. 


In the present study, the compound parabolic collector is assumed 
to be located in Xi’an (34.27N, 108.95E) on June 12, and mounted on 
a horizontal east-west axis and oriented with its aperture plane 
sloping at an angle of 45°. The concentration ratio of the collector is 5, 
the width of its absorber plane is 0.06 m and its length is 2 m. The 
heat collected at the absorber surface is transferred to water flowing 
through two tubes having an outer diameter of 0.018 m and an inner 
diameter of 0.014 m attached to the bottom side. In addition, overall 
loss coefficient U\ 0 is 10.5 W m -2 K -1 , heat-transfer coefficient on 
inside of absorber tube is 230 W m~ 2 K -1 , transmissivity of the cover 


Table 1 

Main assumptions for the CCHP system. 


Environment temperature (°C) 20 

Environment pressure (MPa) 0.10135 

Auxiliary heater outlet temperature (°C) 140 

Heater outlet temperature (°C) 70 

Refrigeration temperature (°C) -5.0 

Turbine isentropic efficiency (%) 85 

Pump isentropic efficiency (%) 70 

Pinch point temperature difference (°C) 10.0 

Approach point temperature difference (°C) 5.0 
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The GA uses fitness function to evaluate adaptability of indi¬ 
vidual without external information in the evolution search. The 
adaptability is expressed by the fitness value. A bigger fitness value 
means a better adaptability subjected to constraints and a better 
viability of the individual. Fitness function which is not constrained 
by definition domain, continuity and differentiability, requires that 
the objective function is defined as the form of non-negative 
maximum. In this optimization, the exergy efficiency is selected as 
the fitness function. 

The GA operators include selection operator, crossover operator 
and mutation operator. Selection operator is responsible for 
selecting the parents to create the next generation of solutions. The 
parent is chosen with a probability based on its fitness. The higher 
the fitness, the higher is the probability of selection. The rank-based 
model is selected for this optimization in the CCHP system. 

Crossover operator is the basic operator for producing new 
chromosomes. It produces new individuals that have some parts of 
both parent’s genetic material. The simple arithmetic crossover is 
applied to this optimization problem due to very simple operation, 
which is presented as follows: 

/ C! = a/i + (1 - a)f 2 ,rr) 

l C 2 = af 2 + (1 - a)/i 1 J 

where a is a random number between 0 and 1,/j and /2 are parents 
individuals which are selected to crossover each other, c\ and c 2 are 
children individuals which are produced by crossover. 

Mutation is needed because even if selection and crossover 
together search new solutions, they tend to cause rapid conver¬ 
gence and there is the danger of loosing potentially useful genetic 
material. The role of mutation in GA has been that of restoring lost 
or unexplored genetic material into the population to prevent the 
premature convergence of the GA to suboptimal solutions. Random 
mutation is adopted to optimize the parameters for the CCHP 
system. It is achieved by selecting individuals from the range of the 
parameter according to mutation probability. 

The steps of GA are made as follows: 

(1) Initialize the population size, crossover probability, mutation 
probability, stop generation, and generate the initial population 
randomly. 

(2) Calculate the fitness of each individual for parent generation, 
and order the fitness. 

(3) Select the individuals from parent generation, and create the 
children generation using crossover and mutation operators. 


Table 2 

Results of simulation for the CCHP system. 



Table 3 

The performance of the CCHP system. 


Turbine work (kW) 

12.730 

Pump 1 work (kW) 

0.079 

Pump 2 work (kW) 

0.327 

Net power output (kW) 

12.324 

Cooling output (kW) 

3.694 

Heating output (kW) 

1.281 

Useful heat gain from solar collector (kW) 

110.708 

Extraction mass flow ratio 

0.5 

Mass flow branching ratio 

0.7 

Heating exergy (kW) 

0.224 

Cooling exergy (kW) 

0.344 

Exergy input from solar collector (kW) 

32.405 

Thermal efficiency (%) 

15.63 

Exergy efficiency (%) 

39.78 


(4) Calculate the fitness of each individual for children generation. 
If the maximum fitness of children generation is less than that 
of parent generation, substitute the maximum fitness of parent 
generation for that of children generation. 

(5) If generation reaches the stop generation, stop the optimiza¬ 
tion. Otherwise go back to step 3. 


4. Results 

The simulation of the CCHP system driven by solar energy was 
carried out using a simulation program written by authors with 
Fortran. Iterative relative convergence error tolerance was 0.01%. 
The thermodynamics properties of the working fluid were calcu¬ 
lated by REFPROP 6.01 [21] developed by the National Institute of 
Standards and Technology of the United States. 

Table 2 shows the thermodynamic state of the CCHP system 
driven by solar energy, and Table 3 shows the results of thermo¬ 
dynamic simulation. 

Fig. 4 shows the variations of the turbine power, heating output, 
cooling output and exergy efficiency with the slope angle of the 
aperture plane for the solar collectors at 10 a.m. (co = -30°J on June 
12 while the other parameters are kept constant as those in Table 1. It 
can be seen that the turbine power, heating output, cooling output 
and exergy efficiency increase first to maximum and then decrease 
as the slope angle of the aperture plane for the solar collectors 



Fig. 4. Effect of the slope angle for the aperture plane on the turbine power, heating 
output, cooling output and exergy efficiency. 
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Fig. 5. Effect of the hour angle on the turbine power, heating output, cooling output 
and exergy efficiency. 

increases. The exergy efficiency reach the maximum when the 
slope angle of the aperture plane is 60°, because the solar collectors 
can absorb the solar radiation effectively at = 60°. 

Fig. 5 shows the variations of the turbine power, heating output, 
cooling output and exergy efficiency with the hour angle at 
f] = 45°on June 12 while the other parameters are kept constant as 
those in Table 1. It is obvious that the turbine power, heating output, 
cooling output and exergy efficiency increase first to maximum and 
then decrease as the hour increases. The maximum exergy effi¬ 
ciency takes place at w = 0°, namely at midday. It is easy to 
understand that the solar radiation is the strongest at midday. 

For the system optimization, the hour angle and the slope angle 
of the aperture plane for the solar collectors are kept to be the 
optimal value of 0°and 60°, and the other parameters of the solar 
collector subsystem are assumed to be constant as those in Table 1. 
Thus, the parameters chosen for optimizing this CCHP system are 
turbine inlet pressure, turbine inlet temperature, extraction pres¬ 
sure, extraction mass flow ratio and mass flow branching ratio 
which could affect the performance of the CCHP system markedly. 
Table 4 shows the condition of the parameter optimization. Table 5 
lists the optimization results for the CCHP system. 

In addition, it is shown that the major output of the system is 
different at any given exergy efficiency. Table 3 shows that at an 
exergy efficiency of 39.78%, 71.24% of the output is net power. Table 


Table 4 

The condition of the parameter optimization. 


Table 5 

The optimization results of the parameters for the CCHP system. 


Turbine inlet pressure (MPa) 

1.183 

Turbine inlet temperature (°C) 

129.62 

Extraction pressure (MPa) 

0.395 

Extraction mass flow ratio 

0.735 

Mass flow branching ratio 

0.231 

Turbine work (kW) 

18.136 

Pump 1 work (kW) 

0.104 

Pump 2 work (kW) 

0.650 

Net power output (kW) 

17.382 

Cooling output (kW) 

4.437 

Heating output (kW) 

77.510 

Heating exergy (kW) 

13.588 

Cooling exergy (kW) 

0.414 

Exergy input from solar collector (kW) 

52.021 

Exergy efficiency (%) 

60.33 


5 shows that at the maximum exergy efficiency of 60.33%, 78.03% of 
the output is heating. 


5. Discussion 

The proposed CCHP system is designed to utilize the solar 
energy to perform the trigeneration based on the crisis of world 
energy in the future. It can supply the electricity, refrigeration and 
heating output simultaneously in terms of different energy 
requirements for an institution (company, plant, hospital, hotel, or 
business center). The amount of the cooling and heating output 
could be regulated by controlling the extraction mass flow ratio and 
mass flow branching ratio to meet the cooling and heating 
requirements of the users in the summer or winter, and the surplus 
electricity can be sold to power grid. In the present study, the 
thermodynamic simulation is conducted to show that the proposed 
CCHP system is feasible theoretically. However, owing to lack of 
experimental data for this CCHP system, the experimental study 
will be carried out to validate the feasibility in the future. In addi¬ 
tion, because of using large numbers of compound parabolic 
collectors to absorb the solar radiation, the investment of the CCHP 
system may be considerable, but this cannot deny the value of the 
CCHP system to solve the energy problem. The economic analysis 
will be conducted to evaluate the feasibility commercially for the 
CCHP system in the coming research. 


6. Conclusion 

In the present study, a new combined cooling, heating and power 
system is proposed. This system is driven by solar energy to reduce 
the consumption of fossil fuel. The system combines Rankine cycle 
and ejector refrigeration cycle, which could produce cooling output, 
heating output and power output simultaneously. The effects of 
hour angle and the slope angle of the aperture plane for the solar 
collectors on the system performance are examined. And the 
parameter optimization is achieved using exergy efficiency as the 
objective function by means of genetic algorithm. 

It is shown that the optimal slope angle of the aperture plane for 
the solar collectors is 60 °at 2 0 a.m. on June 22, and the CCHP system 
can reach the optimal performance with the slope angle of 45° for the 
aperture plane at midday. It is found that the CCHP system has 
a maximum exergy efficiency of 60.33% under the condition of the 
optimal slope angle and hour angle, when turbine inlet pressure, 
turbine inlet temperature, extraction pressure, extraction mass 
flow ratio and mass flow branching ratio are 1.183 MPa, 129.62 °C, 
0.395 MPa, 0.735 and 0.231, respectively. 


Environment temperature (°C) 

20 

Environment pressure (MPa) 

0.10135 

Slope angle for the aperture plane, (°) 

60 

Hour angle, (°) 

0 

Auxiliary heater outlet temperature (°C) 

140 

Heater outlet temperature (°C) 

70 

Refrigeration temperature (°C) 

-5.0 

Turbine isentropic efficiency (%) 

85 

Pump isentropic efficiency (%) 

70 

Pinch point temperature difference (°C) 

10.0 

Approach point temperature difference (°C) 

5.0 

Population size 

50 

Crossover probability 

0.95 

Mutation probability 

0.05 

Stop generation 

200 
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